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Abstract. This paper deals with the Lagrangian analogue of symplectic or point reduc- 
tion by stages. We develop Routh reduction as a reduction technique that preserves the 
Lagrangian nature of the dynamics. To do so we heavily rely on the relation between Routh 
reduction and cotangent symplectic reduction. The main results in this paper are: (i) we 
develop a class of so called magnetic Lagrangian systems and this class has the property 
that it is closed under Routh reduction; (ii) we construct a transformation relating the mag- 
' netic Lagrangian system obtained after two subsequent Routh reductions and the magnetic 

, Lagrangian system obtained after Routh reduction w.r.t. to the full symmetry group. 
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1 Introduction 

It is not a surprise that the bigger part of the recent literature on the geometric description of 
mechanical systems deals in some way or another with symmetry. The benefits of exploiting 
symmetry are indeed self-evident: The equations of motion of a mechanical system which exhibits 
a symmetry can be reduced to a new set of equations with fewer unknowns, possibly easier to 
solve. In particular, the Hamiltonian description of such systems has attracted most of the 
attention in the literature, and the important role played by Poisson manifolds and symplectic 
structures has been strongly emphasized (see e.g. [12] and references therein). Less well-known is 
the process of symmetry reduction for Lagrangian systems. Much like on the Hamiltonian side, 
there are in fact two different paths that lead to two different Lagrangian reduction theories. 
Roughly speaking, the invariance of the Lagrangian leads via Noether's theorem to a set of 
conserved quantities (the momenta). Whether or not one takes these conserved quantities into 
account in the reduction process leads to either the Routh or the Lagrange-Poincare reduction 
method (see e.g. [HElEKiniEKIS]). 

In this paper we deal with Routh's reduction procedure. In a way, one restricts the attention 
to only those solutions of the system with a prescribed value of momentum. The price one 
has to pay is that the new symmetry group (after restriction) is in general only a subgroup of 
the symmetry group of the original Lagrangian system. In [lOj it was pointed out that Routh 
reduction can be interpreted as a special case of symplectic reduction, which moreover preserves 
the Lagrangian nature of the system. 

In the current paper we will investigate those aspects of the theory that are related to re- 
duction in multiple stages. In case the symmetry group of the system has a normal subgroup, 
one can indeed first perform a Routh reduction by means of this subgroup. It is then natural 
to ask whether the reduced system is invariant under some residual, as yet to be determined 
group, and whether a subsequent second Routh reduction leads to results which are equivalent 
to direct reduction by the full group. The answer to the same question, but for the above 
mentioned Lagrange-Poincare reduction theory, is by now well-known [4]. It can be understood 
in the following sense: After each reduction step one remains in the category of so-called 'La- 
grangian systems on Lie algebroids' and symmetry reduction can be phrased in terms of Lie 
algebroid morphisms (see e.g. [5] and the references therein). Equivalently, one may use a con- 
nection to decompose the Lie algebroid structure at each stage. The corresponding category 
is then the so-called category of 'Lagrange-Poincare bundles' from [1], while the corresponding 
'Lagrange-Poincare morphisms' do the reduction. 

It is natural to ask whether or not such a category of systems exists also for Routh reduction 
(it is mentioned as an open problem in [E]). The candidate we propose (in Section [3|) is the class 
of what we call 'magnetic Lagrangian systems'. We show in Sections H] and [5] that this class has 
the property that it is closed under Routh reduction, i.e. that after each step of a reduction in 
stages the reduced system remains in the class of magnetic Lagrangian systems. As in [10], we 
heavily rely on a generalized version of cotangent symplectic reduction. In Section [6] we apply 
the framework of magnetic Lagrangian systems to reduction in several stages, where the first 
stage consists of Routh reduction by a normal subgroup of the overall symmetry group. We end 
the paper with some illustrative examples in the final section. 

2 Preliminaries on Routh reduction 

Definition 1. A Lagrangian system is a pair (Q, L) with Q a manifold and L a function on TQ. 
The manifold Q is the configuration space and L is called the Lagrangian. The dimension of Q 
is called the number of degrees of freedom of the Lagrangian system. A Lagrangian system is 
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of mechanical type if for arbitrary Vq G TqQ, the Lagrangian can be written as 

-^K) = \ {{Vq,Vq))Q-V{q), 

with {{•,■)) Q a Riemannian metric on Q and V a function on Q, called the potential energy of 
the Lagrangian. The function ^ {{vg,Vq))Q is called the kinetic energy. 

We are interested in the Euler-Lagrange equations associated to a Lagrangian system. These 
equations are necessary conditions for a curve q{t) : / C M — )■ Q to extremize the action integral 
Jj L(q(t) , q{t))dt. In a local coordinate chart (g*) on Q the Lagrangian L is a function of (q^,v^) 
and the Euler-Lagrange equations are 

d (dL\ dL ^ i . ^ ^ 

Routh reduction is a reduction technique for Lagrangian systems that are invariant under the 
action of a symmetry group. In classical textbooks, it is a step-by-step procedure that describes 
the construction of a new Lagrangian function with fewer degrees of freedom. 

Before we can formulate Routh reduction, we fix notations for the action of a Lie group on 
an arbitrary manifold M and related concepts. 

Actions of Lie groups and principal bundles. When a Lie group G is given, g denotes 
the Lie algebra of G and exp the exponential map from q to G. The adjoint action of G on g is 
denoted by Ad and the coadjoint action of G on g* by Ad*. 

Definition 2. 

1. A right action of a group G on a manifold M is denoted by '^^'^ : G x M — >■ M; {g, m) i— >■ 

= mg. Throughout this paper we only consider free and proper actions of Lie 
groups on manifolds. 

2. The infinitesimal action is given by : q — t- TmM;£^ ^ d/(ie|o(?7i exp e,^). 

3. For any element ^ S one can consider the fundamental vector field £,m-, defined pointwise 
as iM{m) = ^Am (0- 

4. The map defined pointwise as the dual to ijj^ is denoted by [tp^)* : T*M — >■ g*. 

5. The hfted action ^'^^^ on TM of ^-^^ is given by : G x TM TM;{g,Vm) ^ 

6. The lifted action ^'^**^ on T*M is given by ^'^**^ : G x T*M T*M;{g,am) ^ 

With these notations, ^tm denotes a fundamental vector field on TM determined by the 
lifted action. It follows that ^tm is the complete lift of S^m- 

Every right action ^'^^ gives rise to a left action ^g\m) = ^'^{i(m). We only consider 
right actions. This is not a true restriction since one may reformulate the main results for left 
actions if needed. 

The orbit space M/G of a free and proper action is a manifold and vr : M ^ M/G carries 
the structure of a principal G-bundle. Throughout the paper [mjc denotes a point in the orbit 
space M/G, i.e. it is the orbit through m G M. Every tangent vector in the kernel of Tvr is of 
the form (,M{m) for some € g. These vectors are called vertical and form a distribution which 
we call the vertical distribution Vtt = kerTvr. 

Definition 3. A principal connection on a principal G-bundle is a g-valued 1-form A satisfying 
two conditions: 
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1) it is equivariant, i.e. = Ad^-i • A for any g G, and 

2) for ^ G g arbitrary, A{Cm) = 

The kernel of A determines a right invariant distribution on M which is a complement of 
the vertical distribution. It is therefore called the horizontal distribution of A and is typically 
denoted by H"^ C TM. On the other hand, any right invariant distribution H satisfying 
H © Vtt = TQ determines a principal connection, see [8]. 

We denote by {f,A), where / is a g*-valued function on M, the 1-form on M pointwise 
defined by 

Vm ^ {fim),A{m){u^)) e R. 

The cotangent vector {f,A){m) £ T^M is often denoted by {f{m),A{m)). In particular, if 
/ = // G 0* is constant, then (/i. A) is a 1-form on M. 
Routh reduction. 

Definition 4. Let {Q,L) denote a Lagrangian system and assume that the configuration space 
is equipped with an action . The Lagrangian system {Q,L) is G-invariant if L is invariant 
under the lifted action vj/^'^, i.e. L{"^J'^ (vg)) = L{vq) for arbitrary Vq G TQ and g £ G. 

If a mechanical Lagrangian system is G-invariant then the kinetic energy metric {{■,■)) q and 
the potential energy V are both invariant under the pull-back of for arbitrary g. 

Definition 5. 

1. For a Lagrangian system {Q,L), the Legendre transform ¥L : TQ — )■ T*Q is the fibre 
derivative of L, i.e. for arbitrary Vg,Wg G TQ 



{YL{vg),Wg) = 



L{Vg + UWq). 

M=0 



The Lagrangian is hyperregular if FL is a diffeomorphism. 

2. For a G-invariant Lagrangian system (Q,L), the momentum map Jl ■ TQ — ?> q* is the 
map {il^^Y o FL, i.e. for arbitrary Vq G TQ and ^ G 

{JL{vq),0 = {^L{v,),iQ{q)). 

3. Given a G-invariant mechanical Lagrangian system (Q, L), and a point q in Q. The inertia 
tensor I, is a metric on q defined by = {{iQ{q),'r]Q{<l))) q- 

The momentum map associated to a G-invariant Lagrangian system satisfies the following 
two important properties: 

1. It is conserved along the solutions of the Euler-Lagrange equations, i.e. if q{t) is a solution 
to the Euler-Lagrange equations, then ^ {Jiiqit))) = 0. 

2. It is equivariant w.r.t. the action ^^'^ and the coadjoint action Ad* on g*, i.e. 

Ji(M/J«K)) = Ad;iMvq)) 

for arbitrary g G and Vq G TQ. 
Definition 6. Given a G-invariant Lagrangian system (Q, L) and an arbitrary vector Vq G TQ. 

1. The map Jl\v^ : q ^ q* is the map ^ Jiivq + ^qiq))- 

2. The Lagrangian system {Q, L) is G-regular if JL\vq is a diffeomorphism for every Vq G TQ. 
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Every G-invariant mechanical Lagrangian system is G-regular. To show this, remark that for 
a mechanical Lagrangian system 



Hence JilvqiO = "^i/kq (0) is an affine map, whose linear part is determined by the inertia 



Proposition 1. Consider a G-invariant and G-regular Lagrangian system {Q,L) and fix a re- 
gular value /.i € Q* of the momentum map Jl. Let denote the isotropy subgroup of ^ w.r.t. 



1. The submanifold i^ : (fi) TQ is G ^-invariant and the restricted action of G^ 
on J£^{fi) is free and proper. 

2. The quotient manifold J£^ (i^t)/Gn is diffeomorphic to the fibred product T{Q/G) ^q/g 
QjG^ overQ/G. 

The first statement is a direct consequence of the equivariance of Jl- We postpone a proof 
of the second statement (Proposition [71 page [T5]) . 

We are now ready to describe a preliminary version of Routh reduction where the symmetry 
group G is Abelian and the bundle Q — > Q/G is trivial, so that we can choose a connection A 
with vanishing curvature. Since G is Abelian, we have that the isotropy subgroup G^, for every 
z/ G 0* is the entire group G, and in particular T{Q/G) Xq/g Q/Gu is nothing but T{Q/G). 

Theorem 1 (Routh reduction - Abelian version). Let G be Abelian and let A be a connection 
on Q with vanishing curvature. The Routh reduction procedure of a G-invariant and G-regular 
Lagrangian system (Q, L) consists of the following steps. 

1. Fix a regular value fx of the momentum map Jl and consider the submanifold J£^{fi). 

2. Compute the restriction of the G-invariant function L— {fJ.,A) to the level set Jj^^{iJ,). Let 
L denote its projection to the quotient manifold T{Q/G) = J£^ (fi) / G . 

Every solution of the Euler-Lagrange equations of the Lagrangian system {Q/G,L) is the pro- 
jection of a solution of the Euler-Lagrange equations of the original system (Q, L) with mo- 
mentum /i. Conversely, every solution of the Euler-Lagrange equations of the Lagrangian sys- 
tem (Q, L) with momentum fj, projects to a solution of the Euler-Lagrange equations of the 
system {Q/G, L). 

One may also find a description of global Abelian Routh reduction in [3]. 

Example: the spring pendulum. The system consists of a point particle with mass m 
moving in a horizontal plane, and attached to the origin by means of a spring with spring 
constant k. We choose polar coordinates (r, 6) for this system, so that the mechanical Lagrangian 
is given by L = ^m(r^ + r^^^) — \kr'^. This system is clearly invariant under translations in the 
^-direction, given by ^^(r, 0) = {r,9 -\- a). The momentum map for this action is Jl = mr'^9 
and the Euler-Lagrange equations are 

mr — mrO"^ + A;r = 0, 



metric on g. 



the coadjoint action of G on g*, i.e. g G G^ iff Ad* /j, = /x. 




dt 



Let us fix a regular value 7^ /i = mr'^6 for the momentum map and let A = dO he the standard 
connection with vanishing curvature. The Routhian is obtained from 
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The Euler-Lagrange equation for L is mr 



kr + and solutions of this equations are in 



correspondence to solutions of the Euler-Lagrange equations for L with momentum /i. 

Routh reduction can be extended to more general contexts, including non-Abelian group 
actions |14j . Lagrangians of non-mechanical type [21 [6], Lagrangians invariant up to a total time 
derivative jlOj and Lagrangians that are not G-regular [11]. Additional complications arise in 
these cases: for instance, if the bundle Q ^ Q/G is not trivial, the reduced system is typically 
subjected to an additional force term associated to the curvature of the chosen connection. On 
the other hand, if the symmetry group is not Abelian, the quotient space J^^(/u)/G^ no longer 
has the structure of a tangent bundle; instead it is diffeomorphic to T{Q/G) Xq/g Q/G^. The 
interpretation of the reduced system as a Lagrangian system then requires additional definitions, 
which are postponed to the following section. 

Our main interest in this paper is Routh reduction by stages, where the assumption is that the 
Lagrangian system obtained after applying a first Routh reduction carries additional symmetry, 
so that we can reapply Routh reduction. The system obtained after one Routh reduction is 
formulated on a fibred product of the type T{Q/G) ^q/g Q/G^. This is the fibred product of 
a bundle Q/G^ Q/G with the tangent bundle to the base space. This observation is our main 
motivation to extend Routh reduction to Lagrangian systems defined on such fibred products: 
in order to develop Routh reduction by stages, we have to be able to reduce Lagrangian systems 
obtained after a first reduction. 

Routh reduction is closely related to symplectic reduction on the associated cotangent bundles 
(see [101 I14j). In the next section we will introduce the concept of a "magnetic Lagrangian 
system" and we will emphasize its symplectic formulation. The concept of magnetic Lagrangian 
systems is the analogue of Hamiltonian systems one encounters in magnetic cotangent bundle 
reduction [12]. 

3 Magnetic Lagrangian systems 

A magnetic Lagrangian system is a Lagrangian system with configuration space the total space 
of a bundle e : P — ?• Q ^'^^ where the Lagrangian is independent of the velocities tangent to the 
fibres of e. Additionally the system is subjected to a force term that is of magnetic type. It might 
help to keep in mind that in the case of a Routh reduced Lagrangian system, P corresponds 
to Q/G^ and the fibration e is given by the projection Q/G^ Q/G. 

Definition 7. A magnetic Lagrangian system is a triple (e : P ^ Q,L, B) where e : P — )■ Q is 
a fibre bundle, L is a smooth function on the fibred product TQ xq P and B is a, closed 2- form 
on P. We say that P is the configuration manifold of the system and that L is the Lagrangian. 

A coordinate chart i = l,...,n = dimQ, a = 1,...,/e = dimP — dimQ, adapted 

to the fibration e . P ^ Q determines a coordinate chart {q'^,v^,p^) on TQ Xq P, and the 
Lagrangian L is then a function depending on {q^,v^,p'^). By definition, L is independent of 
the velocities in the fibre coordinates p'^ and therefore it determines a singular Lagrangian 
when interpreted as a function on TP. Locally, the Euler-Lagrange equations for this singular 
Lagrangian are 




i = 1, . . . ,n 



dL 

dp'' 



b 



a = 1, . . . ,k. 



Here we used the following coordinate expression of the 2- form B is ^Bijdq^ A dq^ + Biadq^ A 
dp"" + ^Babdp"" A dp^. These Euler-Lagrange equations have a geometric interpretation (see 
Proposition [2]). First we introduce additional notations. 
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Definition 8. Assume a magnetic Lagrangian system {e : P Q,L,B) is given. 

1. TpQ denotes the fibred product TQ P and {vq,p), where Vq G TQ and p G P such 
that e{p) = g, is a point in TpQ. Similarly, TpQ denotes the fibred product T*Q xq P 
and {aq,p), with e(p) = q, is an arbitrary element in TpQ. 

2. Ve denotes the distribution on P of tangent vectors vertical to e. 

3. e : TP TpQ is the projection fibred over P that maps Vp € TP onto (Te{vp),p) G TpQ. 

4. : TpQ — 7> TQ is the projection that maps {vq,p) G TpQ onto G TQ. 

5. p2 : — > P is the projection that maps {vq,p) G TpQ onto p ^ P. 

6. ei : TpQ ^ T*Q is the projection that maps {aq,p) G TpQ onto aq G T*(5. 

7. £2 : TpQ — )■ P is the projection that maps {aq,p) G TpQ onto p ^ P. 

8. The Legendre transform FL : TpQ TpQ maps {vq,p) G TpQ to {aq,p) G TpQ where 
Uq (zT*Q is determined from 



L(z;g + uwq,p), 

u=0 



for arbitrary i(;g G TgQ. 

9. The energy is a function on TpQ defined by Ei{vq,p) = {¥L{vq,p), {vq,p)) — L{vq,p). 
(Here the contraction of an element {aq,p) G TpQ with {vq,p) G TpQ is defined naturally 

as {{aq,p),{Vq,p)) = {aq,Vq).) 

10. By means of the Legendre transform we can pull-back the presymplectic 2-form e^ajQ + e^^B 
on TpQ to a presymplectic 2-form FL*(eJwQ + eg-B) on TpQ. The latter is denoted by 
qL,b (|jjgj.g = (i^Q, with the Poincare-Cartan 1-form on T*Q and a presymplectic 
2-form is understood to be a closed 2-form, not necessarily of constant rank.) 

Proposition 2. Given a curve p{t) in P, and let j{t) denote the curve in TpQ equal to 
{q{t),p{t)) G TpQ with q{t) = e(p(t)). The curve p{t) in P is a solution to the Euler-Lagrange 
equations for the magnetic Lagrangian system (e : P — ?• Q,L,B) iff ^{t) is a solution to the 
presymplectic equation 

i^(,)J7^'^(7(i)) = -dEMt)). 

Locally, the presymplectic equation coincides with the previously mentioned Euler-Lagrange 
equations, since 

Si''" = d {j^A A dq' + \Btjdq' A dq> + Bi,M A dp' + ^B^dp'' A dp\ 

dE, = „M (14) + - dL = ( fl) - f,dq' - 

\ov^ J ov^ \ov'' J oq^ op"' 

Definition 9. A magnetic Lagrangian system (e : P Q, L, B) 

1) is hyperregular if FL is a diffeomorphism and if the restriction of B to Ve is nondegenerate, 

2) is of mechanical type if L{vq,p) = ^ {{{vq,p), {vq,p)))^_^ — V{p) with ((•, •))^^ is a metric on 
the vector bundle pi : TpQ P and V is a function on P. 

Note that if B is nondegenerate, then the typical fibre of P necessarily has to be even- 
dimensional. In a local coordinate chart the nondegeneracy condition on B is expressed by 
det Bab ^0. 
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Proposition 3. If a magnetic Lagrangian system (e : P Q, L, B) is hyperregular, the 2-form 
^L,B _ ]F^*(^g*(^Q -)- £2-6) determines a symplectic structure on TpQ. 

Proof. Assume that the magnetic Lagrangian system is hyperregular. The nondegeneracy 
of r?-^'^ is easily checked if we work in a coordinate chart adapted to the fibration: 

^""^ = d A dq' + ^B,,dq' A dq^ + Biadq' A dp'' + ^Babdp" A dp\ 

Since ¥L is a diffeomorphism, the 1-forms d[^), dq^ and dp" provide pointwise a basis 
on T*{TpQ). One can now use standard arguments to prove the nondegeneracy. ■ 

We conclude that a hyperregular magnetic Lagrangian system has a symplectic structure 
although the Lagrangian itself is singular when interpreted as a function on TP. The energy El 
is the Hamiltonian. This is an important observation. 

Remark 1. Throughout the paper we only consider magnetic Lagrangian systems that are 
hyperregular. This is not a true restriction. The results remain valid for general magnetic 
Lagrangian systems: instead of relating Routh reduction to symplectic reduction, it is possible 
to relate it to presymplectic reduction [7]. 

Remark 2. We conclude this section with a remark on equivalent magnetic Lagrangian systems. 
Roughly speaking, two Lagrangian systems are equivalent if the resulting dynamics coincide, 
more specifically if they produce the same Euler-Lagrange equations. It is well known that 
for a Lagrangian system the Euler-Lagrange equations do not change when the Lagrangian is 
augmented with a total time derivative of a function on the configuration space. For magnetic 
Lagrangian systems this gauge freedom can be extended. 

Consider a 1-form a along e, i.e. a section of 62 : TpQ P, or in other words a linear 
function on TpQ given by {vg,p) 1— )• {a{p),{vq,p)). When a : P — )• TpQ is composed with 
T*e o ei : TpQ — > T*P, it determines a 1-form on P which, with a slight abuse of notation, 
we denote by e*a. Now consider the function L'{vq,p) = L{vg,p) — {a{p), {vg,p)) and the 
gyroscopic 2-form B' = B + de*a. Together they define a new magnetic Lagrangian system 
{e : P ^ Q, L', B') which is equivalent to (e : P — > Q, L, B). 

Lemma 1. The magnetic Lagrangian systems {e,L',B') and {e,L,B) are equivalent. 

This is a reformulation of a well-known result in classical mechanics, saying that an exact 
gyroscopic force can be taken into account by means of a velocity dependent potential [19j . 

Proof. We show that any solutions to the Euler-Lagrange equations of (e, L', B') is a solution to 
the Euler-Lagrange equations of {e,L,B). We work in a local coordinate neighborhood {q^^p") 
as before. We fix a coordinate expression for a = ai{q^ ,p'')dq^ . The Euler-Lagrange equations 
for (e, L', B') equal, with cf = and p" = v'^ 

d (dL'\ dL' w w a 



dL' 

Qpa 



-^ = -B',y + BUv', a = l,...,k. 



It now remains to substitute the definition of L' and B', i.e. L' = L — aiV^ and 



n' K I dai\ , dai , 

Bij = Bij + [-^ - J , Bi^ = Bia--^, B^i, = Bab- 
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4 Magnetic Lagrangian systems with symmetry 
4.1 Symplectic reduction: a brief introduction 

An overall reference for this section is e.g. [12]. Let (M, w) be a symplectic manifold on which 
G acts on the right, ^^'"'^ : M x G — ^ M. Given a function / : M — ^ g*, then for ^ G g denotes 
the real valued function on M defined by f^{m) = (/(m),^). 

Definition 10. 

1. The action ^''^^ is said to be canonical if (^'^^)*a; = lo for all g € G. 

2. A map J : M ^ g* is a momentum map if i^njU) = —dJ^, for ^ € g arbitrary. 

3. If M is connected, the non-equivariance 1-cocycle a of the momentum map J equals 

a: G^Q*: g ^ J{mg-^) - Ad;., {J {m)), 
where m is arbitrary in M. 

The definition of a is independent of m (see [1]). Recall that a 1-cocycle with values in g* 
statisfies, for g,h G arbitrary, 

a{gh) =CT(g) + Ad*_ia(/i). 

Definition 11. Let o" be a 1-cocycle with values in g*. The affine action of G on g* with 
1-cocycle a is given by {g, fi) ^ Ad*^ -|- a{g~^) for arbitrary £ Q* and g £ G. 

The momentum map J is equivariant with respect to the affine action with 1-cocyle from 
Definition [TUl J{mg) = Ad* J(m) -|- a{g~^). For a fixed element G g*, the Lie group G^ < G 
denotes the isotropy subgroup of ^ w.r.t. the affine action, i.e. g G G^ if fi = Ad*/i -|- a{g~^) or 
equivalently if 

/i - J{mg) = Ad*(^ - J(m)). 

Theorem 2 (symplectic reduction SR). Let {M,uj) be a symplectic manifold, with G acting 
canonically on M. Let J be a momentum map for this action with non-equivariance cocycle a. 
Assume that ^ is a regular value of J, and denote by G^ the isotropy group of ^ under the affine 
action of G on g*. Then the pair {M^,uj^) is a symplectic manifold, with = J~^{iJ,)/G^ 
and with u)^ a 2-form on uniquely determined from z*a; = ttJ^uj^, with i^ : J~^{fi) M and 
vr^: J-H/i)^M^ = J-i(/i)/G^. 

Any Hamiltonian h on M which is invariant under the action of G induces a function 
on satisfying vr*/i^j = The Hamiltonian vector field Xh is tangent to J~^{fi) and the 

corresponding vector field on J~^{fj,) is iTfj^-related to the Hamiltonian vector field X^^ on M^. 



4.2 Invariant magnetic Lagrangian systems and momentum maps 

In order to define invariant magnetic Lagrangian systems, we start from a free and proper action 
of a Lie group G on both P and Q such that they commute with e: for arbitrary g £ G, 

6o^; = xI/«oe. (1) 

In other words, G acts on e : P — )• Q by bundle automorphisms. The projections of the principal 
bundles are denoted by ir^ : Q Q/G and vr^ : P — >■ P/G. These actions induce lifted right 
actions on TpQ and TpQ: 
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1) ^-f := {T^'^{vg),^^{p)) for {vg,p) € TpQ and geG, 

2) ^^^^{ag,p) := (r*M/J_,(a,),M/f (p)) for {a„p) G T*pQ and 5 e G. 
Definition 12. A magnetic Lagrangian system (e : P — > Q,L,B) is G-invariant if 

1) B is invariant under the action of {"^^YB = B for all G G, 

2) the Lagrangian L is invariant under the lifted action ^i'^p^ on TpQ, i.e. L o ^I/J^*^ = L for 
all geG. 

Let (e : P ^ Q,L,B)hea G-invariant hyperregular magnetic Lagrangian system and consider 
the symplectic 2-form e^LOQ + on TpQ. 

Proposition 4. The action ^I'^p^ is canonical and the Legendre transform ¥L : TpQ — > TpQ 

is equivariant, i.e. FL o ^J^*^ = ^J^*^ o FL. As a consequence the action "^^p^ is canonical 
w.r.t. $7^'^. 

Proof. Note that ei o ^p*^ = T*'^^_^ o ei and 62 o ^rp*^ = o £2. Since the action = 
T*^'^_i on T*Q is canonical w.r.t. uq and since S is a G-invariant 2-form on P, we conclude that 
\^TpQ jg canonical. To show the equivariance of FL we use its definition: let g ^ G, Vq, Wq G TQ 
be arbitrary, then 

(FL (^p«(i;„p)) , M'f = L(^p'3(i;, + uwq)) 



d_ 

du 



L{Vq + UWq) = {¥L{Vq,p), {Wq,p)). 

u=0 



In order to define a momentum map, we introduce a concept that is similar to the so-called 
Sg-potential in [T2] . 

Definition 13. A 0*-valued function 5 on P is a Sg-potential for the magnetic 2-form B if 
is^pB = d5^ for any ^ G g. 

Recall that 5^ is defined pointwise as S^{p) = {d{p),(,). If B is G-invariant and admits 
a ;Bg-potential 6, then it follows that d{{"i?g)*d^ — (^Adg^) = 0. If P is connected, this induces 
a g*-valued 1-cocycle on G: 

as: G^Q*; g^ - Ad;_, • 5{p), 

with p arbitrary in P. The infinitesimal version of this 1-cocycle is denoted by 'Ssi^.^'n) — 
— {Teas{^,),r]). It is a 2-cocycle on the Lie algebra and satisfies £5(^,77) = —B{^p,r]p) — S^^^^j = 
^p{6r)) — For the following proposition, recall from Definition [2] that {ip'^)* : T*Q g* is 

the dual of the infinitesimal action of G on Q. 

Proposition 5. The map Js = (ip^)* oei — 6oe2 is a momentum map for the symplectic manifold 
{TpQ,e\ujQ + €2^) with non- equivariance 1-cocycle —as- Due to the equivariance of¥L, the 
map Jl^s = Js o FL is a momentum map for the symplectic manifold {TpQ, 17^'^) with non- 
equivariance 1-cocycle — (T5. 

Proof. The map : T*Q — ?> g* is an equivariant momentum map for the symplectic 

manifold {T*Q,uiq). It is straightforward that the combined map Js = o ei — 5 o £2 is 

a momentum map for the lifted action on TpQ w.r.t. the symplectic form e^wg + e^^B with 
non-equivariance cocycle — (T5. The rest of the statement follows by construction. ■ 

We conclude that the symplectic structures {TpQ,^}'^'^) and {TpQ,elujQ -|- eg^) associated 
to a G-invariant and hyperregular magnetic Lagrangian system with a jBg-potential 5 admit 
a momentum map with cocycle —as and are amenable to symplectic reduction. 
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5 Routh reduction for magnetic Lagrangian systems 
5.1 Magnetic cotangent bundle reduction 

In this section we study the reduction of the symplectic manifold {TpQ, e\ojQ + e^B) determined 
from the magnetic 2-form ;B of a G-invariant hyperregular magnetic Lagrangian system with 
i3g-potential 5. Prom Proposition [5l it follows that the function Js is a momentum map with 
non-equivariance cocycle —a^. As usual, denotes the isotropy subgroup of /.t G g* for the 
affine action. 

Proposition 6. Fix a connection A'^ on the bundle tt^ : Q — )• Q/G and a regular value fi of 
the momentum map Js- 

1. g(^G^lff^l + 5{pg) = Ml{^l + 5{p)). 

2. = €*A^ is a principal connection on vr^ : P P/G. 

3. The 1-form {fi + 6,A^) on P is G ^-invariant. 

4. The 2-form B + + 5,A^)) is G ^-invariant and reducible to a 2-form on P/G^. 

5. The quotient manifold J^^(fi)/Gf^ is diffeomorphic to Tp^^ {Q/G). 

Proof. Recall the pointwise definition of {ijl + 5, A^): given Vp € TP, then + 5, A^){p){vp) := 
(/i + 5{p),A^{p){vp)). Below, {n + 5{p),A^{p)) G T*P denotes the cotangent vector (/i + 
5,A^){p). We continue with the proof. 

1. This is straightforward from the definition of cr^. 

2. The 0*-valued 1-form A^ determines a principal connection if it satisfies A^{£,p) = ^ for 
ah e G 0, and if (^f = Ad^-i • A^ for all 5 G G: 

A''{ip){p) = AQ{e{p)){Te{ip{p))) = A'^ieipmQieip))) = ^ and 
{^PyA^ = {^'g)*e*A^ = e^'^fYA^ = Ad^-i • e*A^ = Ad^-i • A^ . 

3. The pull-back 1-form {^g)*{^l + b.A^) equals 

{^i + 5{pg),{{^PrAP){p)) = (Ad*(/. + 5(p)),Ad,-i -^^(p)). 

4. That B -\- d{{ix + S^A^)) is G^-invariant is a straightforward consequence of 3. It is 
projectable to P/G^ since the contraction with any fundamental vector field with ^ G 
vanishes: 

i^p {d{{fi + 6, A^)) +B)= C(.p{{^ + 6, A^)) -d{{^ + 6,i)) + dS^ 

= (Cp(Ai + S),-^^) + + 5, ad_5 • A^) - d6^ + dS^ 
= {^p{6) - adlfi - adp, A^) = 0, 

where we used the fact that ^ G g^, or equivalently ad|/i = i^T^s- 

5. Similar to cotangent bundle reduction, we use a 'shift map' to construct the required 
diffeomorphism. Recall that Vtt'^ denotes the bundle of tangent vectors vertical to vr'^ : Q — >■ 
Q/G. The subbundle V^tt^ of T*Q is defined as the annihilator of Vtt^. Below we introduce 
a shift map (p"^ between and VpTT^ = V^tt'^ ^qP and we show that it is equivariant w.r.t. 
to actions of obtained by restriction of ^^^'3 J~i(^) and VpTr'^ . This is sufficient for (p'^ 
to project to a diffeomorphism [c/)'^] between the quotient spaces. This provides us with the 
desired diffeomorphism because the quotient of VpTT*^ is well known: VpTr^ /G^ = Tp,^ {Q/G). 
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Definition of the shift map. In the following we use a slight abuse of notations: if we 
write (/i + S{p), (q)) for some p G P and q = e{p), then this is the cotangent vector in T*Q 
determined by Vq t-^ {fi + 6{p), A'^{q){vq)). Let 0^ be the map 



J5 ^ ^P^Q'^ iaq,p) ^ {ag- {ii + 5{p),A^{q)),p). 
It is well defined, i.e. (t>'^{ctq,p) G VpTT^ for any {aq,p) in the level set J^^{pL) since 

(<K,P),(Cg(9),P)) = K,eQ(g)) - {lJi + 6{p),A'^{q){iQ{q))) 
= {Jsiaq,p),O-{f^,O=0. 

Equivariance of the shift map. For arbitrary g E G^, 

= (r*^^-iK) - {f^ + d{pg),A'^{qg)),pg) 

V^^_, (a, - {f, + S{pg),Ad^-i-A'^{q))) ,pg) 

(a, - {f^ + 6{p),A'^{q))) ,pg) = v&f^ (</>^(a„p)) . 



The projection of (f>-^ is thus well defined, and it is denoted by [</>^], i.e. [0^] : Jg ^{^)/G^ — >■ 

Definition 14. S is the 2-form on P/G^ obtained after reducing B + d{{ii + 5,A^)). 

In the following proposition we put a symplectic structure on the fibred product Tpi^ (Q/G) 
= T*{Q/G) Xq/g P/G^. The notations are similar to those on TpQ. 

Definition 15. We define the projections 

1) 6i:T*/^^(Q/G)^r*(Q/G), 

2) 62:r*/^^(Q/G)^P/G^. 

The canonical symplectic 2-form on T*{Q/G) is denoted by ojq/g- 

Completely analogous to the construction on TpQ, we introduce e\ujQjQ + e^^B as a 2-form 

onr*/^^(g/G). 



M, = J^\^)/G, 



T*{Q/G) 



Figure 1. Cotangent bundle reduetion (notations as in Theorem [2] with M ~ TpQ). 
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Theorem 3 (Generalized magnetic cotangent bundle reduction). Let fj, denote a regular value of 
the momentum map Js for the symplectic structure (TpQ, e^(^Q+e2i3) with isotropy subgroup Gfj_, 
and let he a principal connection on Q — )• Q/G. Then the symplectically reduced manifold 
{M^,uj^) is symplectomorphic to {Tp^^ {Q/G),elujQ/Q + e^^B)), with symplectomorphism [0^]. 

Although this result is a straightforward extension of standard magnetic cotangent bundle 
reduction in |12j . we include a proof because its method turns out to be useful later on. 

Proof. Let iq denote the injection V^tt'^ TpQ and 7r° the projection V^tt^ Tp^Q^{Q/G). 
We know that [4>'^] is a diffeomorphism, and therefore it only remains to show that the 2-form 
^i^Q/G + is pull-backed to under this map (and therefore, implicitly elujq/Q + will 
be nondegenerate) . 

We use the fact that is uniquely determined by i*^ujQ = 7r*w^, with : J^^ifJ.) — )■ TpQ 
the natural inclusion and vr^ : J^^ifJ.) the projection to the quotient space. Due to the 

uniqueness property, it is therefore sufficient to show that 

<.([0^]*(ela;Q/G + ~elB)) = i* (e^c^g + e^B). (2) 

The left-hand side of ([2]). Since [<^jf]o7r^ = 7rJ^°0^ the left-hand side of the above equation 
is {^loct>p*[e\u:Q/G + ~elB). 

The right-hand side of ([2]). We first show the equality of two 1-forms: 

i^ielOQ) and {^i)* (il{e\eQ + el{ii + 6, A""))). 

Let {aq,p) G J^^{fi), V(^aq,p) ^ TiTpQ) a tangent vector to J^^{fi) and Vq £ TgQ denotes the 
projection of V(^ag,p)- Then, the first 1-form equals 

^^l(ei^Q)(ag:P)(^K,p)) = ^Q(ag)(^ei(^K,p))) = ("g,^?)- 
The second 1-form evaluated on this tangent vector gives 

(<)*(i5(el^Q + 4{l^ + 5,^^)))(a„p)(V(,,,p)) 

= OQiaq - (M + 5(p),^«(g)))(r6i(r(</.^)(V(,^,rt))) + {^L + S{p),A<^{q){vq)) 

The two 1-forms are identical and after taking the exterior derivative and adding £2^, we obtain 
that the right-hand side of ([2]) equals 

i^ieluQ + e*2B) = {^^y{i*o{elu;Q + e^(d(M + 5,^^) + B))). 

Since (i) (€2 o TfO)*^ = (£2 o ioYidifi + 5,^^) + B), and (ii) {h o vrO)*^^/^ = i^elOg, one easily 
verifies that the left-hand side equals the right-hand side in (l2|). ■ 

We conclude with a result on the behavior of a symplectomorphism under symplectic re- 
duction. Assume two symplectic manifolds (M, 0) and (M',f]') and a symplectomorphism 
/ : M — )■ M' (i.e. a map for which f*il.' = il.) are given. We assume in addition that both 
M and M' are equipped with a canonical free and proper action of G. Let J : M ^ q* and 
J' : M' — 7> Q* denote corresponding momentum maps for these actions on M and M' respec- 
tively. We say that / is equivariant if f{mg) = f{m)g for arbitrary m (z P, g (z G. Note that the 
non-equivariance cocycles for J and J' are equal up to a coboundary. Without loss of generality 
we may assume that f*J' = J and that the non-equivariance cocycles coincide. This in turn 
guarantees that the affine actions on q* coincide and that the isotropy groups of an element 
fJ- £ Q* coincide for both affine actions. Finally, fix a regular value /i € 0* for both J and J' . 
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Theorem 4. If f is an equivariant symplectic diffeomorphism M — t- M' , such that J' = J o f , 
then under symplectic reduction, the symplectic manifolds {Mfj_,0,^) and are symplec- 

tically diffeomorphic under the map 

[f,\: M,^M'^; Hg, ^ [fim)]G,. 

Proof. This is a straightforward result. Since / is a diffeomorphism for which J' = J o /, 
the restriction of / to J~^{ii) determines a diffeomorphism from J~^{ij) to J'~^(/u). The 
equivariance imphes that reduces to a diffeomorphism [/^] from = J~^(/i)/G^ to M'^ = 
J'~^{^)/G^. It is our purpose to show that [ffj\*^'^ = ^f_i, or since both vr^ and vr^ are projec- 
tions, that vr*f]^ = f^{TT'*Cl'^). The determining property for Q,^ is vr*f]^ = (and similarly 
for ^2^). From diagram chasing we have that = /^(i^J^')- Then 

Tr;n^ = i;n = f;{i';n') = f;w;n'^) = 

since tt^ o = [f^] o vr^^ by definition. This concludes the proof. ■ 
5.2 Routh reduction for magnetic Lagrangian systems 

In Proposition [5] we have introduced the momentum map Jj^ g = o FL for a G-invariant 
hyperregular magnetic Lagrangian system (e : P — )• Q,L,B), i.e. 

{JlA^„p),0 = {¥L{vg,p),{^Q{q),p)) -5^{p). 

We know from Proposition [J] that the Legendre transform FL is equivariant. After restriction 
to J£^g{fi) it reduces to a symplectic diffeomorphism between the symplectically reduced spaces 
J^^(/i)/G^ and J^^{^)/G^ (see also Theorem The following diagram summarizes these 
previous observations. 

{TpQ, n^'^) {T*pQ, eluQ + e^B) 



SR 



SR+[(/)' 



.^1 



['A;f]o[(FL)M] 

{Jl^{^l)/G^,^l^) ^ iT*p/^JQ/G),ilu;Q/G + -e*2B) 

Figure 2. Diagram relating tangent and cotangent reduction (SR = symplectic reduction). 

We are now only three steps away from a description of Routh reduction for magnetic La- 
grangian systems. The final goal is to describe the symplectic reduced manifold (J^^(^)/G^, 0^) 
as a symplectic manifold associated with a new magnetic Lagrangian system, with Lagrangian 
say Lfj_. 

1. The first step is the construction of a diffeomorphism : Tp/G^{Q/G) — )■ J^^(/i)/G^. 

This is crucial because the Hamiltonian dynamics determined by the Euler-Lagrange equations 
reduce to Hamiltonian dynamics on a manifold of the form Tp/G (Q/G). If we can characterize 
this reduced dynamics as the Euler-Lagrange equations of a magnetic Lagrangian system with 
configuration manifold -P/G^ — >■ Q/G, we have developed a Routh reduction technique for mag- 
netic Lagrangian systems. This characterization consists of the two remaining steps mentioned 
before: 
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2. We show that the composition [4>'^\ o [(FL)^] o equals the fibre derivative of some 
Lagrangian function L on Tp/Q^{Q/G), and this imphes that the symplectic 2-form A*r2^ 
on Tp/Q^{Q/G) is of the form 

FZ*(e>Q/G + eie), 

which is the symplectic structure associated to the reduced magnetic Lagrangian system 
{P/G^^Q/G,L,By, 

3. We show that the reduction of the energy El in the symplectic reduction scheme is pre- 
cisely the energy Ej^ on Tp/Q^^{Q/G). This guarantees that the Euler-Lagrange equations 
associated to the reduced Lagrangian L are related to the Euler-Lagrange equations for L. 

The new reduced Lagrangian L is what we call the Routhian (it is often alternatively denoted 
by R or R^). This is summarized in the diagram in Fig. [3l 

r n 



SR+A;1 



{Tp/G^Q/G),n^'^) {T* {Q/G),ilu;Q/G + ~e*2B) 



Figure 3. A symplectic view on Routh reduction. 



Step 1: the definition of : Tp/Q^{Q/G) — >■ Jj^\{^)/G^. First consider the map 
: J^l{fi)/G^ — )• Tp/Q^{Q/G) which is obtained from the G^-invariance of the map 



J£}{fi) ^ Tp/G^iQ/G); {v„p) ^ {T7TQ{vg), \p] 



G n 



We will define the map as the inverse of 11^. In general however, 11^ is not invertible. The 
following definition is the analogue of G-regularity for a classical Lagrangian. 

Definition 16. The Lagrangian L of a G-invariant magnetic Lagrangian system is called G- 
regular if the map JL,&\(vq,p) ^ ~^ 0*;C JL,s{vq + iQ{<l)-,P) is a diffeomorphism for all 

iVqiP) G TpQ. 

Every mechanical G-invariant magnetic Lagrangian system is G-regular because JL,s\(vq,p) is 
an affine map modeled on the metric on g obtained from the kinetic energy metric. 

Proposition 7. 11^ is a diffeomorphism if the Lagrangian is G-regular. 

Proof. We construct the inverse for 11^. Choose an element (^[gjg, Mg^) in TpjG^iQ /G) and 
fix a point {vq,p) £ TpQ that projects onto [pIg^)) i-e. Tir'^^Vq) = V[g] and p G [p]g^- Due 

to the G-regularity of L, there exists a unique element ^n g such that JL,s{vq + Cq{q),p) = Z^- 
The orbit of the point {vq + S,Q{q),p) under the G^-action determines an element in J^J(/i)/G^. 
We will show that this construction is independent of the chosen point {vq,p), and therefore 
defines a map A^j : TpjG^{Q/G) — > J£\{fi)/Gfj_ which is the inverse to 11^. 

To show that the construction is independent of the chosen point {vq,p) in TpQ, we choose 
any other point {Vq,p') in TpQ projecting onto [p]g^)- The point is of the form (v'q,p') = 

{"^g {vq + 'nQ{Q))iP9) for some g € G^ and G 0. We now repeat the previous construction 
applied to {v'q,p'): we consider a unique element ^' G g for which Ji^si^'g^ {vq + flQio)) + 
i'Q{l9)iP9) — ^- Due to the equivariance of Jl.s-, we conclude that 

// = JL,5{^1^{vq + riQ{q)) + i'Q{qg),pg) = Ad* • Ji^sivq + (r? + Mgi')Q{q),p) - asig"^). 
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Since g € G^, we may conclude that JL,5{vq + {r]+ Adg^')Q{q) , p) = fi and therefore r]+Adg^' = ^. 

This imphes that [Vg + £,q{q)iP') = '^li^^i'^q + ^q{q)}P) foi' 9 ^ G'^ such that p' = pg. This 
concludes the proof. ■ 

The previous proposition guarantees that for a given G-regular Lagrangian and a bundle 
adapted coordinate chart {x^,g^) in Q — )• Q/G and {x'^,g^,p°') in P — > Q, the functions 
{x^ , , g"^ , p^) determine a coordinate chart in J^^(/i). This is useful in the following. 

Step 2: the reduced symplectomorphism. We will now introduce a new Lagrangian L 
on Tp/Q^{Q/G) with the property that its Legendre transform ¥L : Tp/Q^{Q/G) — )■ Tp^^ (Q/G) 

is precisely the reduced Legendre transform [(FL)^], more specifically ¥L = [(P^] o [(FL)^] o A^. 

As is already described in |10j . the structure of the Lagrangian L is completely determined by 
the cotangent bundle reduction scheme, and in particular by the shift map (f)-^ from J^^(^) — > 
Vpir'^ . To make this statement more precise, we remark that [4>'^] o [(FL)^] is the quotient map 
oicp-^oFL. Thus, if [{vg,p)]G, is arbitrary in J^j(/i)/G^, then [cP^] o [{¥L)M^g,p)]G,) is the 
quotient of 

¥L{v„p)-{{fi + 6{p),A'^{q)),p). 

The latter is the fibre derivative of a new function on TpQ, namely L{vq, p) — {iJ,+5{p) , {q){vq)) . 
This function is clearly G^-invariant. 

Definition 17. The function L is defined as the pull-back under of the quotient map of the 
restriction to J^^ifJ') of the function 

^ L{Vq,p) - {fi + 6{p),A^{q){Vq)) . 

Lemma 2. ¥L = [(f>-^] o [(FL)^] o A^,. 

Proof. Fix elements (fjg]^, [pjc^) G Tp/Q^{Q/G) and fix a representative {vq,p) € Ji\{fJ-) in 
the orbit A^(z;[g]^, [p]g^) £ /^^i-- definition of the maps involved, we have 

([0;^] o [(FL)^] o A^) [p]gJ = (vr° o cl>^){¥L{vq,p)) 

= iTl{¥L{vq,p) - {{^1 + 5{p),A'^{q)),p)). 

Next we study the fibre derivative of the map L. Fix a point (""^[gj^, [p]g^) find we compute 

d 



(^L{V[q\G^ WgJ, (w^Mg, NgJ/ 



M=0 



We construct a curve n i-)- C('w) in projects onto the curve u i— A^(f[^]jj + 

^^[(jIg) [p]g^) in J^^(;u)/G^ such that C(O) = {vq,p) and C(0) is vertical to the projection 
P2 ° ifj. ■ JL\iP') ~^ ^ (recall that p2 : TpQ ^ P is the projection onto the second factor in 
the fibred product). 

The existence of such a curve is best shown in a specific coordinate chart (see above) for 
JJ^^{p): {x^ , x'^ , g'^ , p"") where the index r = 1, . . . , dim Q/G, m = l,...,dimG and a = 
l,...,dimP — dimQ. Note that {x'^,x'^) is a coordinate chart in T[Q/G). In these coordi- 
nates, we let {vq,p) = (x'q, Xq, 5™,pg) and w^g^^ = {xq,Wq). We define the curve C{u) to be the 
curve u i-> {xq,Xq + uw^^g^ ,Pq). Then the tangent vector to Tpi{C,{Q)) is the vertical lift of 
some Wq € TqQ with Tir'^lwq) = w^gj^. 
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Finally, from the definition of L and the fibre derivative FL we obtain 



F^(%]g, WgJ, (t^Mc WgJ) = ^ 



M=0 



since ¥L{vq,p) — {{fi + S{p),A'^{q)),p) € VpTr'^. This concludes the proof. ■ 

Step 3: the reduced energy Hamiltonian. The third and last step concerns the specific 
reduced dynamics. We have to relate the energy El on TpQ to the energy of the Routhian L 
on Tp/Q^{Q/G). In the following Lemma we again use the notations from Theorem [2] applied 
to the reduction of the symplectic structure on {TpQ, Q^''^). 

Lemma 3. The energy is the reduced Hamiltonian, i.e. it satisfies: 

iU^o7r,yE-^ = i;EL, 
with ovr^ : J^j(/x) J^j(/i)/G^ Tp/G^{Q/G) and i^ : J^jlpt) ^ TpQ. 
Proof. Let {vq,p) G J^^{^), such that (H^ o 'n^){vq,p) = (^'[gjQ, [v\g^)- Then 

i*^EL{Vq,p) = {¥L{Vq,p), {Vq,p)) - L{Vq,p) 

= {{^■^o{¥L)^){Vq,p),{Vq,p)) - {L{Vq,p) - {fl + 6 (p) , A'^ (q) (Vq))) 

= {{[cp-^] o [(FL)^] o A^)(t-[,]^, [p]g,), (%]^, [ph,)) - HgJ. 
Using the result from Lemma [2] this concludes the proof. ■ 

Routh reduction. The previous three steps are summarized in the following theorem. 

Theorem 5 (Routh reduction for magnetic Lagrangian systems). Let (e : P ^ Q,L, B) be a hy- 
perregular, G-invariant and G-regular magnetic Lagrangian system and let 5 he a Bg-potential 
of the magnetic term B with 1-cocycle as- 

1. Let fj, Q* be a regular value of the momentum map J^^g and let G^ be the isotropy 
subgroup of fi w.r.t. the affine action on g* with 1-cocycle —(Js, i-c. g € if and only if 
II = Ad*g fi - as{g~^). 

2. Fix a connection A^ on vr*^ : Q — )> Q/G and let A^ be the corresponding connection 
on P ^ P/G. Compute the restriction of the G^-invariant function L{vq,p) — (/x + 
6{p), A^ {q){vq)) to J£^g{fi) and let L be its quotient to Tp/Q^{Q/G) = J^J(/i)/G^. 

3. Compute B as the projection to P/G^ of the 2-form B + d{fi + 6, A^). 

4. Consider the magnetic Lagrangian system: (e^ : P/G^ — )• Q/G,L,B). 

This reduced magnetic Lagrangian system is hyperregular and every solution p{t) ^ P to the 
Euler-Lagrange equations for (e, L, B) with momentum fj, projects under P — > P/G^ to a solution 
of the Euler-Lagrange equations for (e^ : PjG^ — >■ Q/G,L,B). Conversely, every solution to 
the Euler-Lagrange equations for (e^ : PjG^ — > Q /G, L,B) is the projection of a solution to the 
Euler-Lagrange equations for (e, L, B) with momentum fi. 

It is possible to say more about the structure of B and its relation to the connection A'^. 
This and reconstruction aspects fall out of the scope of this paper. We refer to [11] where 
these topics are described in more detail. They carry over to this more general framework in 
a straightforward way. 
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5.3 Reduction of magnetic Lagrangian systems on Lie groups 

Consider a magnetic Lagrangian system on P = Q = G, i.e. the configuration space is a Lie 
group G. We start from a function ^ on g and with it we associate a Lagrangian L on TG 
by right multiphcation L{g,Vg) = £{vgg~^). By definition L is invariant under the right action 
of G on itself. We assume that a magnetic 2-form B is given which is invariant under right 
multiphcation and admits a fig-potential (5 : G — > g*. We first rephrase some definitions in this 
specific setting. 

1. The 1-cocycle cj^ : G — t- g* satisfies (Ts{g) = 6{hg~^) — Ad* ^i6{h) , for arbitrary h ^ G. If we 
let h = e, then as{g) = S{g~^) — Ad*-i6{e) or equivalently, Ad*as{g) = Ad*g6{g^^) — 6{e). 
Similarly if we let h = g, then as{g) = 6{e) — Ad*-i6{g). Since 6 is determined up to 
a constant, we may assume without loss of generality that (5(e) = 0. 

2. The associated 2-cocycle T,s{(,,r]) = (.ci^ri) — ^[^,r)]- 

3. We use the right identification of TG with Gxg, i.e. {g,Vg) is mapped to {g^Vgg"^) G Gxg. 
The right action of G on TG equals right multiplication in the first factor of G x g under 
this identification. 

4. We use the Maurer-Cartan principal connection on G G/G: A{g){vg) = g~^Vg. In the 
right identification, the connection corresponds to the map (5, ^) € G x g i— > Ad^-i^ € g. 

5. The momentum map Jl^s ■ TG g* equals JL,5{g,^9) = ^<i*g¥i{^)-5{g), for (5, ^) G Gxg 
arbitrary. If L is G-regular then is invertible, i.e. there exists a function X • 0* ~^ 
such that ¥£{xiiy)) = v. 

6. The affine action on g* is (5^, /u) 1-^ Ad*/x — as{g~'^)- 

7. The isotropy group G^ consists of group elements g such that + 6{g) = Ad*/i. 

8. The quotient G/G^ (right coset space) can be identified with O^, i.e. [gjc^ £ G/Gf^ is 
mapped onto u = Ad*„i/i — crs{g) = Ad*_i(/i + S{g)). A tangent vector to G/G^ at [g]Gf^ 
which is the projection of S^g is mapped to a vector = — ad^z^ + i^T,s in g*. 

Lemma 4. The 2-form B + d{{fi + 5, A)) reduces to the Kirillov-Kostant-Souriau symplectic 
2-form on = G/G^. 

Proof. Let g be arbitrary, and let Vg = ^g, Wg = rjg be two tangent vectors in TgG with ?/ G g 
arbitrary. Note that Vg = S,Q{g) with ^' = Ad^-i^, and similarly Wg = 'r]Q{g) with r]' = Adg-irj. 
Then 

B{g){vg,wg) = B{g){i'G{g),v'G{g)) = Vcik'M- 
On the other hand 

d((/. + <5, ^))(5)(z;„ = + -5, ^'))(5) - r/^((/U + 5, 0)(5) - (/^ + 5(5), [e', r/']). 

Before continuing, we compute the equivariance of the 2-cocycle Tis{S,' 

d d 
^5i^',v') = -^{crsig^^ e^P s^g),rj') = --j-{as{g~^ exp s^) + Ad*^p_^^gas{g),v') 

= -^{as{g-^) + Ad*gas{exp s^) + AdlAd*^^_^^(7s{g), Adg-i7]) 
Summarizing, we have 

{B + diif, + 6, A))) {g){vg,wg) = + v')){9) - + le, V]) 
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This 2-form is reducible to a 2-form on . If we use the isomorphism G /Gfj_ — t- introduced 
above, the 2-form reduces to B{i'){i', v') = (z>, rf) with r/ € such that v' = — ad*z^ + ij^S^. ■ 

We conclude by computing the Routhian L as a function on O^. By definition, it equals the 
reduction to of L — {fi + 6, A)\j^ g=^. In the right identification, the level set Jl,s{9: '^g) = A* 
is precisely Ad*_i(/i + 5{g)) = F^(^), with = Vg. If we set v = Ad*_i(/i + 5{g)) € Of,, the 
fixed momentum condition is = xi^)- The Routhian L(i') becomes in the right identification 

By application of the chain rule, it easily follows that {dL[u),i'') = — {i'' , xi'^)) ■ The reduced 
Euler-Lagrange equations are 

iyB{v) = dL{v) or z> = -ad*(j,)i/ + (3) 

For later purpose, we remark that for a left action and £ originating from a left invariant 
Lagrangian, the reduced equations are z> = a.(l*,-.u — i^(^y-)Tjs (here the 1-cocycle satisfies as{g) = 



6 Routh reduction by stages 

In reduction by stages, we study the reduction of a G- invariant system (symplectic or La- 
grangian) under the action of the full group G and under the induced action w.r.t. a normal 
subgroup K <\ G. We shall adopt as far as possible the notations used in [12]. A detailed 
construction of the following definitions is found in this reference. 

Definition 18. 

1. The Lie-algebra of -ftr is and i denotes the injection i : q with dual i* : q* ^ ^* . 

2. The group G acts on ^ by restriction of the adjoint action. The induced action of G on 
^* is denoted by the same symbol Ad* : G x — > J^* . 

3. /X denotes an element in g* and u G ^* . Then is the isotropy subgroup of /i under the 
Ad*-action of G on g*; Gi, is the isotropy subgroup of i' under the Ad*-action of G on ^* 
obtained as the dual of the restricted Ad-action of G on K; and K^, is the isotropy of i' 
w.r.t. to standard coadjoint action of K on A*. These groups satisfy Gi, (iK = Ky and Ky 
is normal in Gy. 

4. Qi, and denote the Lie algebras of Gy and Ky respectively. Gy denotes the quotient 
group Gy/Ky and its Lie algebra equals Qy = Qy/^y 

5. The projections onto the quotient groups are denoted by r : G — > G = G/K and Vy : Gy ^ 
Gy, and on the level of the Lie algebra: r' : q ^ q = g/^ and r'y : Qy ^ q,^. The inclusion 
map Gy ^ G induces a map ky : Qy ^ g, with its dual /c* : g* — )• g*. 

6. p denotes an element in g*. 

In [12] symplectic reduction by stages is performed under the condition of a so-called 'stages 
hypothesis'. An element € g* is said to satisfy the stages hypothesis if for any fi' € g* 
satisfying = /i'|t^ = v and /i|g,^ = iJ,'\g^ = u, there exists an element k G Ky and g G {Gy)y 
such that Ad^^/i' = The stages hypothesis is a condition on a chosen momentum value and 
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depends only on the symmetry group G. It was akeady clear in [12] that the hypothesis is 
automatically satisfied if G is a central extension or if G is a semi-direct product group. In the 
recent contribution [16] it has been pointed out that the hypothesis is in fact always satisfied, 
and that it can be taken out of the reduction by stages statements altogether. Taking advantage 
of this result, in this paper, we will not make further reference to the stages hypothesis. 

6.1 Symplectic reduction by stages 

Theorem 6 (Symplectic reduction by stages [E]). Let {M,u;) be a symplectic manifold with 
a canonical G-action with an equivariant momentum map Jq- 

1. Fix a regular value fi (z Q* of the momentum map and perform symplectic reduction to 
obtain the symplectic manifold (M^,;^^). 

2. The restriction of the action to K is canonical and the map Jk = i* ° Jg ■ M ^ K* 
determines an equivariant momentum map for this induced action. Fix a regular value v 
of Jk and perform symplectic reduction to obtain the symplectic manifold {My,ojy). 

3. The level set J^^{v) is G,y -invariant. 

4. The group G^, = G^/K,^ acts on Mp by projecting the restricted action of G^, on Jj^^ii'). 
This induced action ^^^'^ is free, proper and canonical on [My^ijjy). Assume that Ky is 
connected. 

5. Fix an element v in g* such that the restriction ofv\^^ equals v\^^. There is a well-defined 
momentum map Jq^ : My — >■ g* for the induced action '^^•^ . This momentum map is 
determined from Jq and v and has a non-equivariance cocycle: Jq^ satisfies 

{f'uT ° Jg^ oTTy = kl o Jg o iy - i?. 

6. Let p (z qI be a regular value for the momentum map Jq^ and let {Gy)p be the isotropy 
subgroup of p w.r.t. the affine action of Gy on g*. Perform symplectic reduction to obtain 
the symplectic manifold {{My)p, {ujy)p) with {My)p = J^^ {p)/{Gy)p. 

If p is chosen such that (r'y)* p = p\g^ — v, then there exists a symplectic diffeomorphism 

F : (M^,W^) ^ {{My)p,{LOy)p). 

For our purpose it is also important to understand the reduction of a G-invariant Hamil- 
tonian h on M. We assume that all conditions in Theorem [6] are satisfied. First note that, 
by definition of the momentum maps, we have an inclusion jp of Jq^{p) in J^{u). Recall 
that we use vr^ for the projection Jf^^iv) — > My. It was shown in jl2j that the image of 
T^v ° jp '■ Jc^ifJ') ~^ My is contained in Jq^{p). Moreover, this map is equivariant w.r.t. the 
action of Gp on Jq^{p) and {Gy)p on J^^{p) (this makes sense, since Gp projects to a subset 
of {Gy)p). The quotient of T^v°jp is the symplectic diffeomorphism F mentioned in the previous 
theorem (see also Fig. H]). 

Let /i be a G-invariant Hamiltonian on M and let hp be the function on Mp obtained from 
vr* /i^ = i*^i_H. On the other hand we let hy be the function satisfying iTyhy = i*h. This function is 
G,y-invariant: hy{[m]K^[g]K^) = h{iy{mg)) = h{iy{m)g), with m G J^^iv) and g ^ Gy arbitrary. 
Note that (tTj, o jp)*hy = i*^h. 

The Hamiltonian hy is a Gj^-invariant function on {My,ujy). Applying the second symplectic 
reduction to this manifold, we obtain a new reduced Hamiltonian {hy)p on {My)p. 

Proposition 8. F*{{hy)p) = hp. 
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(M,) 



Figure 4. Commuting diagram relating the different reduced symplectic manifolds. 



Proof. We rely on the commuting diagram in Fig. [5| 

7r;{F*{{K),)) = {7r,oj^rK = i;{h). 
This uniquely characterizes F*{{hi,)p) as the function h^. ■ 

6.2 Routh reduction by stages 

Routh reduction by stages is symplectic reduction by stages applied to the symplectic struc- 
ture of the initial Lagrangian system. In this section we show that the symplectic structures 
and energy hamiltonians in the different stages can in fact be associated to specific magnetic 
Lagrangians systems, and eventually gives us Routh reduction by stages. The symplectic re- 
duction by stages then provides us a diffeomorphism relating the solutions of the different 
Euler-Lagrange equations for the Lagrangian systems in the final stages. 

We start with a hyperregular Lagrangian L, invariant under the action of a Lie group G. We 
assume that this Lagrangian satisfies a regularity condition which is more stringent than mere 
G-regularity. 

Definition 19. The Lagrangian L is said to be G-hyperregular if for any Vq S TQ and any 
subspace < g with injection i' : ^' ^ g, the mapping i'* o J^jt,^ o i' : ^' ^ j^'* defined by 
e ^ i'*{JL{vq + {i'{0)QiQ))) is invertible. 

Lagrangians of mechanical type are G-hyperregular. Let K be a normal subgroup of G, K the 
Lie algebra of K and i : q the canonical injection. Due to the hyper regularity the invariant 
Lagrangian L is both G- and ET- hyperregular, and both G- and iiT-invariant. By definition of Jl, 
the map i* o Ji = i* o [ip'^^)* o FL is the momentum map for the if-action. 

Theorem 7 (Routh reduction by stages). Assume {Q,L) is a hyperregular, G -hyperregular and 
G-invariant Lagrangian system. Let K denote a normal subgroup of G. 

1. Let fJ. & Q* be a regular value of the momentum map Jl and a G -connection on Q. Let 
{Q/G^ — )• Q/G, Lq, Bq) be the magnetic Lagrangian system obtained by performing Routh 
reduction with respect to G. 

2. Fix a regular value v ^ ^* of the momentum map i* o for the K-action and a K- 
connection on Q. Assume that A^ is G-equivariant w.r.t. the action of G on Q and ^* . 
Consider the magnetic Lagrangian system {Q/Ky Q/K,Li,Bi) obtained by performing 
Routh reduction with respect to K . 

3. Gy acts on Q/J'Cjy and Q/K by projecting the induced action of Gy on Q. These induced 
actions are free and proper. 

4. Assume that is connected. Fix an element v ^ such that i>\i^^ = z^Ij^^,. Then 
the magnetic Lagrangian system {QjKy Q/ L'C, Li, Bi) is Gy-invariant, Gy-regular and 
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admits a BiQu -potential 5i entirely determined by the choice of y. The potential satisfies, 
for arbitrary q & Q 

Let Ji denote the momentum map associated with 6i. 

5. Fix a regular value p ^ qI, for the momentum map J\ and let {Gy)p be the isotropy sub- 
group of p w.r.t. the affine action of Gy on Qjy. Fix a G^/- connection on Q/K. Consider 
the magnetic Lagrangian system {{Q/Ky)/{Gy)p — > {Q/ K)/Gu), L2, B2) obtained by per- 
forming Routh reduction with respected to Gy. 

If p is chosen such that {r'^)* p = /ijg,^ — v, then every solution -y{t) G Q/G^ to the Euler- 
Lagrange equations for [Q/G^ — ?> Q/G, Lq, Bq) is mapped to a solution in [Q/Ky)/{Gi,)p to the 
Euler-Lagrange equations for {{Q/Ky)/{Gy)p {Q/K)/Gy,L2.,B2)- Conversely, a solution in 
{Q/Ky)/{Gy)p to the Euler-Lagrange equations for {(Q/Ky)/(Gi,)p — )■ {Q/K)/Gu, L2, B2) is the 
projection of a solution in Q/G^ to the Euler-Lagrange equations for [Q/G^ — )■ Q/G,Lq,Bo). 

Proof. 1 and 2 are obtained by applying Routh reduction. 3 follows from |12[ p. 152]: we know 
that the quotient groups G^, = Gy/Ki, acts in a free and proper way on the quotient space QjKy. 
The group Gy is a subgroup of G and acts freely and properly and Q/K. We now show 4. 
G,y-Invariance of the Routh reduced system {Q/Ky — > Q/K, Li,Bi). 

Lemma 5. If the connection is chosen such that it is equivariant w.r.t. the action of the full 
group G, i.e. if 

{^^rA' = Ad^-^A\ 

then the magnetic Lagrangian system (Q/Ky — )■ Q/ K, Li, Bi) is Gy-invariant and 61 is a BiQy- 
potential. 

Proof. We first show that Li is G,y-invariant. For that purpose, we choose an arbitrary g € Gy 
and let g € Gy he a representative. Similar we choose a point (v^qj^, [q]K^) S Tq/x^{Q/ K) such 
that it is the projection of Vq £ {i o Ji)~^(u) c TQ. By definition of the quotient action on 
Tqix^{Q / K), the action of g on an element {v[q]j^, [l]Ky) equals the projection of Vqg. We now 
check the invariance of Li at an arbitrary point in Tq j i^^{Q / K): 

L, [qW)) = - {i^,A\qa){Tm<^{vq))) 

= L{vq) - {Adl^.u, A\q){Vq)) = Li {v^q^^, [q]K^) . 

Next, we check the Gjy-invariance of Bi. Recall that Bi is the projection to Q/Ky of the 2- form 
d{v,A^) on Q. We first consider the equivariance of this 2-form under Gy. Let g € Gy be 
arbitrary, then 

{^^Y{d{u,A^)) = d{Ml^,v,A^) = d{u,A^). 

We thus obtain G,y-invariance for d{v,A^), and we may conclude that {^f^'^ fBi = Bi holds 
on Q/Ky. 

The third and final step is the definition of the ,Si0,y-potential. We consider an element 
? = ^ Qu = Qu/^u and let ^ £ Qy he a representative. Then, by definition of Bi, the 1-form 
^^Q/K projection to Q/Ky of the 1-form i^Qd{{i', A^)) on Q (i.e. projects to £,q/k„)- 

Again we concentrate on the 1-form on Q: 



i^^d{u, A') = {{y, A^)) - d{i^^ {v, A^)) . 
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Since {ij,A^) is Gi,-invariant, we conclude that i^^d{i', A^) = —d(^i^^{i',A^)). The exact 1-from 
on the right gives a strong hint of the structure of the 3,y-potential. Assume now that we fixed 
an element T' £ qI such that = ulj^^. 
The function 5 on Q, defined by 

-6^iq) = {,.,AHq){CQ{q)))-{V,0 

is our candidate for the fiig^-potential. This statement makes sense provided that 5^ projects 
to a function on Q/K^ and that it only depends on the equivalence class ^ = ^ + of ^ € g^. 
The latter is a straightforward consequence of the fact that A^ is a principal iT-connection. The 
-fCjy-invariance is more involved, and we rely on a result in [12]. For any k in K^, we have 

-6^{qk) = {u,A\qk){iQ{qk))) - {V,^ = A\q){{AdkOQ{q))) " (^,0- 

Therefore 6^ is constant on the orbits of K,^ in Q if {i',A^{q){{S, — Adfc^)Q(g))) vanishes for 
all k. To show this we introduce a function / on K,y given by f{k) = {u, A^{q){{^ — AdkC)Q{q))) 
and we use similar arguments as in [121 p. 156]. If we can show that /(e) = 0, df\e = and 
7(^1^2) = f{ki) + f{k2) for arbitrary ki^2 G Ky, we may conclude that / = (since is 
assumed connected). 

The first condition /(e) =0 is trivial. To check the second condition: let k € .^jy be arbitrary, 
then 

dfUK) = {u,A\q){-s.d^OQ{(D)) = -(^,ad,0 = -(ad:,0 = 0. 

Above, we have used the fact that Kjj is normal in and that, as a consequence, the Lie 
bracket [k,^] is in Therefore the contraction of the corresponding fundamental vector field 
with A^ is precisely [k,^]. Next, we check the third condition and compute f{kik2). Given the 
identity 

i - AduMk.i = e - Adfc,^ + Adfc,(^ - Adk,0 
and the fact that ki G K,^, 

fikik2) = {u,A\qm-Adk,k,OQ{q))) 

= W,A\qm - M,,OQiq))) + {iy,A\q)iAd,,{C - Adk^OgiQ))) 
= f{ki) + (Ad^^z., A\qm - Adfc,e)Q('7))) = fih) + f{k2). 

This completes the proof: the 0*-valued function 5 is shown to be projectable to a 0*-valued 
function on QjKy. This is the sought-after potential b\: for arbitrary q ^ Q^we have 

KnSrMK.)) = -{^l^^^ ok,r{{v,A\q))) +V. ■ 

Symplectic structure of {Q/Ky — > Q/ K, Li,Bi) and symplectic reduction by stages. 

Lemma 6. Apply symplectic reduction by stages to the symplectic structure associated to the 
G-invariant Lagrangian system {Q,L). Identify the symplectically reduced manifold Mi, with 
the symplectic structure on Tq/j^^{Q/K) induced by the magnetic Lagrangian system {Q/K^ — > 
Q/K,Li,Bi). Then: 

1. The action xJ/^Q/^fi^ W/^) of G^ on Tq^j^^(Q/ K) is precisely the induced action on the first 
reduced space Mi, in symplectic reduction by stages. 

2. For a chosen u G g*, the momentum map Ji : Tqip^^^{Q/K) g* associated with the mag- 
netic Lagrangian system {Q/Ky — )• Q/K,Li,Bi) corresponds to the induced momentum 
map Jq^ from symplectic reduction by stages. 
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Proof. 1. The momentum map for the K action is precisely Jk '■= Jl, with Jl : TQ g*. 
By definition, the induced action of on J^^{u) / is obtained by projecting the action 
of on J^^(z^). If we take into account that we reahze the quotient manifold Jj^{v)/Ky as 
Tqij^^{Q/K), the induced action on Tqij^^{Q/K) is obtained by projection of the action of Gy 
on TQ under the projection TQ Tqi^^iQ / K). This is precisely the action we have introduced 
above. 

2. The induced momentum map Jq^ is defined in the following way (we consider it directly 
as a function on Tqix^{Q / K) instead of on My): 

{jG,XmKM]K.),~0 = {JL{v,),ky{0) - = (FLK),eQ) - (4) 

where ^ G 0,^ is arbitrary and projects to ^ € Qy, Vq projects to W\K^)- By definition of 

the momentum map of the magnetic Lagrangian system [Q/Ky — t- Q/K, Li,Bi), we have 

{Ji{v[q],, [q]K.),C) = {^M^I,]k^ [q]k.), HOq/k, \p]k..)) - {5i)^{[p]kJ. (5) 

We now show that the right-hand side of ([5]) equals the right-hand side of We therefore use 
the definition of FLi and 6i as being the projection of maps upstairs: 

(FLi(z;[,]^, [q]K^), {{Oq/k, [qU..]) = {¥L{vq),^Q{q)) - {i.,A\q){CQ{q))), 

i5l)X[Q]K.) = SdP) = -W,A\q){^Qiq))) + (17, e). ■ 

Finally, before we can reapply Routh reduction for the second stage, we need to check that 
Li is Gj^-regular. 

Lemma 7. The magnetic Lagrangian system (Q/Ky — ?> Q/ K, Li, Bi) is Gy-regular. 
Proof. We have to show that, for any (t'[g]^, [q\kv) the map 

is invertible. Let Vq determine a tangent vector in J^^^if) representing {v^q^^,[q\K^)- Let rj 
denote an arbitrary element in g*. Due to the assumed G-hyperregularity, there is a unique 
^ ^ Qy such that k*{JL{vq + Cq('?))) = ('^[/)*^ + ^- The projection = rl{^) of defines the 
inverse element for rj, since it is such that 

KiMv, + ^q(.q))) -y = {r'yY Jxiv^q^^ + Iqik) = KTv- ■ 

Symplectic and Routh reduction by stages. Summarizing the above lemmas, we con- 
clude that the magnetic Lagrangian system {Q/Ky — )• Q/K, Li,Bi) is amenable to Routh re- 
duction and that the symplectic structure and momentum map associated to this Lagrangian 
system correspond to the symplectic structure and momentum map encountered in symplec- 
tic reduction by stages. If p is chosen such that the compatibility relation {r'y)*p = plg^ — i> 
holds, then from symplectic reduction by stages we have that the symplectic structures asso- 
ciated to {Q/Gy, Q/G,Lq,Bq) and {{Q / Ky)/{Gy)p {Q/ K)/Gy, L2, B2) are symplectically 
diffeomorphic by means of the symplectic diffeomorphism F introduced earlier. From Proposi- 
tion [8] it follows that F*El^ = Elq and therefore the corresponding Hamiltonian vector fields 
are F-related. We define a map r : Q/G^ {Q/Ky)/{Gy)p as T{[q]G^) = Wq]k„'\(^q ) • The 
map is well-defined since G^ is a subgroup of Gy and since ry{Gp) C {Gy)p. 

Lemma 8. The symplectic diffeomorphism F is fibred over r. 
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Tq/g^Q/G) ^ T^Q/K.)KG.),m/K)/G,) 



QIG^ {Q/K,)/{G,)p 

Figure 5. Fibratioii of tlie symplectic diffeomorpliisni F. 

Proof. We recall the definition of the map F: fix an element (^'[5]^, Mg^) and let Vq G J£^{fi) 
be a representative. The point F{v[q-^^, [g]G^) is obtained by taking the consecutive quotients 
of Vq. In particular, the component of the final quotient in the configuration space {Q/ K^)/{Gu) 
of the magnetic Lagrangian system, is precisely the image r. ■ 

Since -F is a diffeomorphism, r is onto. The Hamiltonian vector field on Tq/q^{Q/G) and 
'^{Q/K^)/{G^)piiQ/-^)/^'^) -F-related. Their integral curves project onto solutions of the 
Euler-Lagrange equations. This concludes the proof of Theorem [71 ■ 



7 Examples 
7.1 Elroy's Beanie 

This system appears in e.g. [13] . It consists of two planar rigid bodies that are connected in their 
center of mass. The system moves in the plane and it is subject to some conservative force with 
potential V. The configuration space is SE{2) x , with coordinates {x,y,6,tl^). Here {x,y) is 
the position of the center of mass, 6 is the rotation of the first rigid body, and -0 the relative 
rotation of the second body w.r.t. the first. The kinetic energy of the system is 5i?(2)-invariant 
and we will suppose that the potential is invariant as well. This means in fact that only the 
relative position of the two bodies matters for the dynamics of the system. The Lagrangian is 
of the form 

L = Kn{±^ + 2/2) + + + - Viij). 

The Euler-Lagrange equations of the system are, in normal form, 

x = o, j/ = o, e = ^v\ ;p = -h±Jiv'. 

The symmetry group and the principal connection. An element of SE{2), the special 
Euclidean group, can be represented by a matrix of the form 

(6) 



cos a 


— smt) 


X 


sin 9 


cos 6 


y 








1 



The identity of the group is (x = 0, y = 0, 9 = 0) and the multiplication is given by 

{xi,yi,ei) * {X2,y2,d2) = (x2 cos 6^1 -1/2 sin 6*1 +xi,X2sin6'i + 7/2 cos 6*1 +yi,6i + 92). 
The matrices 




ei = u u u , 62 = u u i , 63 
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form a basis for the Lie algebra, for which [61,62] = 0, [61,63] = 62 and [62,63] = —61. A corre- 
sponding basis of fundamental vector fields on Q is 

d d Odd 

OX ay ox oy oti 

and a basis of invariant vector fields is 

„ d . ^ d „ ■ n 9 ^ d „ d 

ei = cosfc*— + sm6'— , 62 = - sm t*— + cos 6*— , 63 = —. 

ox oy ox oy off 

One can easily verify that the Lagrangian is invariant under the S'i?(2)-action. There is a trivial 
principal connection on P = Q = SE{2) x 5^ — )• Q/SE{2) = S"^, which locally takes the form 

{dx + yd9)ei + {dy — xd9)e2 + ^^63. 

The momentum map Jl is given by: 

Jl = mxe^ + mije^ + {m{xy — yx) + Ii9 + l2{6 + V'))^^- 

In what follows we perform two Routh reductions on the Lagrangian. The first reduction is 
done w.r.t. the full symmetry group SE{2), and the second reduction w.r.t. the Abelian normal 
subgroup M?. 

Full reduction. Let = fiie^ + /-f26^ + /^36^ be a generic element in g*. An element 
^ = .^^61 + ^^62 + ^'^63 of the isotropy algebra satisfies 

CV2 = 0, eVi = 0, eV2 - eVi = 0. 

So if we suppose that ^1 and /i2 do not both vanish - we will set ^1 = 1 from now on - then 
a typical element of is ^ = ^^(61 + 112^2)- Since g^j is 1-dimensional, is of course Abelian. 
A convenient way to describe the manifold P/G^ = SE{2)/G^ locally is by considering a coor- 
dinate transformation {x',y',6') in the group coordinates such that the vector field associated 
to an element in g^ becomes a coordinate vector field: in the new coordinates, we should have 
dx' = dx + /i2<9y. This is obtained by the following transformation 

x' = x, y' = y - fJ'2X, 6' = 9. 

Then clearly {y',9) is a coordinate chart on SE{2)/G^. And simultaneously, we have that 
{y', 9, ip) is a coordinate chart on the reduced configuration manifold P/G^ = {SE{2) x S^)/G^, 
and the fibration P/G^ — )• Q/G = is locally represented by {y',9,ip) 1-^ {ip). For a more 
systematic treatment on appropriate coordinate changes, we refer to [6]. 

We now compute the Routhian Lq and the 2-form Bq. Following [19], a convenient way 
to compute the (unreduced) Routhian for mechanical Lagrangians is by using 2{Lq + V) = 
( — PxX — Pyij — 9pe + i^Pip) j-'^^^y where pi is the momentum in the ith coordinate. We have: 



2(Lo + V)= + y2) _ j^q2 _ ^ + l2{9 + i;)^^ 



-m{±^ + y2) _ (/^ + /2)^2 ^ j^^2^ 



2 



1 ^ ( - {xfi2 -y)- hip) 

— (1 + M ) - ^ J— J ^ + 

m h + h 

hh _^2 _^ /^3 - {xii2 - y) ^ (/is - {x^i2 - y)f 



h+h h+h h+ h 
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In the last step we have left out some constant terms. The reduced Lagrangian is then ob- 
tained by taking the quotient w.r.t. the action of G^. This is done by applying the coordinate 
transformation introduced above. We get: 

Lo = + h^i^ - ( vii.) + lii^i±yt] , 

° 2/1 +/2 h+h \ ^^ 2 I1+/2 J' 

which is clearly independent of x' . The 2- form Bq is obtained by reducing the 2-form 

d{{dx + ydO) + ^2{dy - xdO) + ^JL^dO) = d{y - ^.2x) A dO. 

Using the coordinate change we get Bq = dy' A dO. 

In this example the Routhian Lq depends on the velocity corresponding to the coordinate ijj 
on (S*^, but is independent of the velocities corresponding to the two remaining coordinates (y', 6). 
With the above, the reduced Euler-Lagrange equations take the form 

y =0, 

1 



h + h 



Note that the second order equation in ip decouples from the first order equations, and that 
these two first order equations are the momentum equations rewritten in normal form. 

Abelian reduction. We now perform first Routh reduction w.r.t. the Abelian symmetry 
group of translations in the x and y direction. Let us denote the symmetry group by iiT = 
and study the quotient spaces. We will use the same notations as before: the Lie algebra elements 
ei, 62 denote a basis for the subalgebra of X in g. The momentum map for this action is now 
Jk = i* o Ji^ = mxe^ + mije^ . We choose v to be the projection of the momentum we had 
used in the full reduction: let v = + /i2e^ € g*. Since K is Abelian, = K and the quotient 
space is SE{2) x S^/M^ = 5^ x S^. If we choose = dxei + dye2 to be the trivial connection, 
we simply get Bi = 0. The Routhian Li can now be obtained from 

2{Li + V) = ( -p^x - pyy + pgO + p^Tp) 

-mx^ - mf + hO"^ + hie + Tpf) = h0^ + h{0 + i^f, 

^ Jl ('^) 

where we ignored again some constant terms. The Routh reduced system is now a standard 
Lagrangian system on x with Lagrangian Li = ^hO"^ + ^h{d + ip)"^ — V{ip) (see the 
paragraphs on Abelian Routh reduction). Its equations of motion are 

ii iih 

For this example there is actually no second stage: the group Gy = GyjKy is the trivial 
one {e}, and the vector space gjy = Qvj^v is only the zero vector. So, there is no second 
momentum map to take into account, and there is no further symmetry to quotient out. 

In the reduction by stages process we have not made use of ^i-^. We now show that the two 
ways of reducing the system are equivalent. 

Equivalence between direct reduction and reduction by stages. Let us compute the 
diffeomorphism F for this example. Here, it is a map J]p^{^)/G^ — ?• J^^(u)/K that is obtained 
by projection of the inclusion map J£ (fi) (v). The latter equals, in coordinates 

{x, y, e, ip, Tp) ^ (x, y, 9,tp,e = -^-^-^(^3 - (2:^2 - v) - hip),ip 
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The map is reducible, and after taking the quotient it becomes 



F: Jl\i,)/G,^J],\u)/K- 



1 



{f^3 + y' - hij),^^ 



h + h 



This diffeomorphism maps the G^-reduced system on the J^jy-reduced system, as is obvious 
from the respective equations of motion. 

7.2 Rigid bodies on the Heisenberg group 

As a second example of Routh reduction by stages, we discuss the dynamics of a rigid body 
immersed in a potential flow with circulation |20j . We assume that the body is circular, and in 
this case the equations of motion are given by 



Px 


= r 


'-Vy 


, where 


Px 


= M 


Vx 


Py. 








Py_ 




Vy_ 



(7) 



Here F represents the circulation and M is a (non-diagonal) mass matrix, which incorporates the 
inertia and added masses of the body. The right-hand side of the equations of motion represents 
the so-called Kutta-Joukowski lift force, a gyroscopic force due to circulation [9l [T7] . 

While this system is extremely easy to integrate, it nevertheless exhibits all the interesting 
geometric characteristics of more complicated examples. As we show below, the configuration 
space for this system is the Heisenberg group, arguably the simplest non-trivial central extension 
group, and the procedure of reduction by stages demonstrated here can be applied equally well 
to more complicated central extensions, such as the oscillator group (describing the dynamics of 
rigid bodies of arbitrary cross section in circulatory flow) and the Bott-Virasoro group describing 
the KdV equation. 

In this context, the Heisenberg group H is the Euclidian space M^, equipped with the multi- 
plication 



using the definition of H given in [18]. The coordinates (x, y) describe the center of the circular 
disc, while the coordinate s corresponds to the flux of the fluid around the body. The conjugate 
momentum corresponding to s will be the circulation of the fluid. 

The Heisenberg group can alternatively be described as the central extension of M? con- 
structed by means of the cocycle B{{x,y), {x' ,y')) = \{xy' — yx'). We now introduce the 
infinitesimal cocycle C : x ^ M, given by 

C{{v.:^,Vy),{v'^,v'y)) =V^v'y-VyV'^. 

For the relation between B and C, see for instance |12] . The Lie algebra \) oi H can then be 
identified with with the bracket 

\{vx,Vy,v),{y'^,v'y,v')\ = {0,0,C{{vx,Vy),{v'^,Vy))) = {0,0,VxVy - Vyv'J. 

The dual Lie algebra [)* can again be identified with M^, with coordinates {px,Py,p) and duality 
pairing {{px,Py,p), {vx,Vy,v)) = PxVx + PyVy + pv. 

We now induce the following quadratic Lagrangian £ on f): 




i{Vx,Vy,v) = - [Vx,Vy\ M 
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and we define L on by left extension: L{g,g) = £{g ^g), or explicitly 



L{x,y,s;x,y,s) = - [x,y]'. 



X 




y. 





(8) 



The Euler-Poincare equations obtained from I are given by 



d 
di 



Px 




-Vy 


Py 


= P 


Vx 


.P. 




_ _ 



After setting p = F, we obtain the equations ([7]). These equations also coincide with the Euler- 
Lagrange equations obtained from the Lagrangian ([8]). 

Reduction with respect to M. The center of H is the normal subgroup isomorphic to M 
which consists of all elements of the form (0, 0, s), where s G M. We first perform Routh reduction 
with respect to the left action of this subgroup on H . On the principal bundle H — )• H/M. = 
we consider the connection one-form given at the identity by A{e){vx.,Vy,v) = v, and extended 
to the whole of H by left translation. Explicitly, we have 



A{x,y,s) 



ds — - {xdy — ydx) 



Since the structure group M is Abelian, the curvature of A is given hy B = dA = —dx A dy. 

Similarly, the momentum map Jl ■ TH R for the R-action on the tangent bundle TH is 
given by JL{x,y, z;x,y, z) = s — {xy — yx)/2, so that consists of all points {x,y, s;x,y, s) 
with s = r + (xy — yx)/2. The isotropy subgroup Rr = R acts on this level set by translations 
in the s-direction, so that the reduced velocity space is given by 

J~i(r)/Rr = TR^ 

The symplectic form on the reduced space can easily be computed, and is explicitly given by 

Adx t\dx + B{dx Ady + dy A dx) + Cdij A dy - Fdx A dy, 

where A, B, C are the entries of the mass matrix M. The last term of the symplectic form, 
—Tdx A dy, is the curvature term of the connection, paired with F € R. Finally, a quick 
computation shows that the reduced Lagrangian is just the kinetic energy Lagrangian on TR^: 



Li{x,y;x,y) = ^ [x,y\M 



(9) 



up to constant terms. 

Second reduction. We now perform reduction with respect to the remaining symmetry 
group, iJ/R = R^, using the results from Section 15.31 We have a left invariant magnetic 
Lagrangian system on the group R^, with Lagrangian ([9|) and magnetic form By = —Tdx A dy. 
The potential 5 : R^ — > R^ corresponding to the latter is given by 



5{x,y) 



-y 



and the momentum map is therefore J2{x,y,x,y) = 
2-cocycle of the momentum map is Y^s = B-p- 

The affine action of R^ on itself is given by {x,y) ■ {px,Py 



{x,yY^ — 5{x,y). The non-equivariance 



{px - Ty,py + Fx). If we fix 



a momentum value (A, n) € R , the isotropy group 



of the affine action consists of just the 



zero element, and consequentially the twice-reduced space J2 ^(A, /i)/R^_^ is nothing but 
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The reduced Euler-Lagrange equations ([3]) in the case of a left action, assume the fohowing 
form 



d 
It 



Px 

Py. 



with {pxiPyY^ = ^{vx^Vy)'^ , and these are nothing but the equations ([7]). 
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